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Abstract 

We analyze the production and subsequent decay of the neutral MSSM Higgs 
bosons {h = h°, H°, A ) mediated by flavor changing neutral currents (FCNC) 
in the LHC collider. We have computed the /i-production cross-section times the 
FCNC branching ratio, a(pp — > h — > qq') = o{pp —>h)x B(h — > qq'), in the LHC 
focusing on the strongly-interacting FCNC sector. Here qq' is an electrically neutral 
pair of quarks of different flavors, the dominant modes being those containing a 
heavy quark: tc or bs. We determine the maximum production rates for each of 
these modes and identify the relevant regions of the MSSM parameter space, after 
taking into account the severe restrictions imposed by low energy FCNC processes. 
The analysis of cr(pp — > h — > qq') singles out regions of the MSSM parameter space 
different from those obtained by maximizing only the branching ratio, due to non- 
trivial correlations between the parameters that maximize/minimize each isolated 
factor. The production rates for the bs channel can be huge for a FCNC process 
(0.1 — lpb), but its detection can be problematic. The production rates for the tc 
channel are more modest (1CP 3 — 10~ 2 pb), but its detection should be easier due 
to the clear-cut top quark signature. A few thousand tc events could be collected 
in the highest luminosity phase of the LHC, with no counterpart in the SM. 



1 Introduction 



The search for physics beyond the Standard Model (SM) is a very relevant, if not the most 
important, endeavor within the big experimental program scheduled for the forthcoming 
Large Hadron Collider (LHC) experiment at CERN [1,2]. There are several favorite 
searching lines on which to concentrate, but undoubtedly the most relevant one (due 
to its central role in most extensions of the SM) is the physics of the Higgs boson(s) 
with all its potential physical manifestations. Supersymmetry (SUSY) [3-6] is certainly 
related to Higgs boson physics, and at the same time it may convey plenty of additional 
phenomenology. Ever since its inception, SUSY has been one of the most cherished 
candidates for physics beyond the SM, and as such it will be scrutinized in great detail 
at the LHC. It is no exaggeration to affirm that the LHC will either prove or disprove 
the existence of SUSY, at least in its most beloved low-energy realization, namely the one 
which is needed to solve the longstanding naturalness problem in the Higgs sector of the 
SM [4]. Indeed, if SUSY is realized around the TeV scale, the LHC experiments shall be 
able to directly produce the SUSY particles for masses smaller than a few TeV [7,8]. On 
the other hand, the presence of SUSY may also be tested indirectly through the quantum 
effects of the supersymmetric particles. For one thing it has been known since long ago 
that SUSY particles may produce large virtual effects on Higgs boson observables 1 . 

Within the general strategy based on detecting indirect effects of the new physics, 
Flavor Changing Neutral Currents (FCNC) play a very special role. In the SM they 
are completely absent at the tree-level [23]. At one loop, however, the FCNC effects are 
possible in the SM, but then the contributions from the new particles enter on equal footing 
with those from SM particles. Therefore, it is not inconceivable that in certain regions 
of the parameter space the new physical effects may well dominate the SM contributions, 
thus providing a unique signature of physics beyond the SM. This is particularly so when 
the SM one-loop effects, even though non-vanishing, turn out to be highly suppressed. In 
such situations the sole observation of these FCNC processes would be instant evidence of 
new physics. A most dramatic example of this kind of appealing scenarios occurs within 
the FCNC physics of the SM Higgs boson (H SM ) interactions with quarks. For example, 
the FCNC vertex Hs M tc may lead (at one loop) to such rare decays as t — > H SM c or 
Hsm — > t c, depending on the mass of Hsm (At the moment both possibilities are still 
open because the LEP bounds amount to a mass for Hsm above 114.4 GeV [24].). Both 
of these modes are extremely suppressed at one loop, with branching ratios of order 10~ 14 
or less, hence 10 orders of magnitude below other more conventional (and relatively well 
measured) FCNC processes like b — > S7 [24]. The rareness of the FCNC Higgs boson decay 
modes is borne out by direct calculation [25], and it can also be understood from physical 
arguments based on dimensional analysis, power counting, CKM matrix elements and 
dynamical features [26-28]. Therefore, these processes are an ideal laboratory to look for 
non-standard interactions superimposed onto the SM ones. Similar considerations apply 
to the FCNC processes associated to the Hbs vertex, but in this case it is more difficult 
to pin down the phenomenological signatures. Some work along these lines has already 

1 See e.g. [9-21] and references therein. For a review see e.g. [22]. 
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been done, both in the MSSM [26, 29-36] and in the general two-Higgs-doublet model 
(2HDM) [27,28,32,37], and also in other extensions of the SM- see [38] for a review. Up 
to now, the main effort has been concentrated in computing the FCNC decay modes at 
one-loop within the new physics, and also in getting a realistic estimate of the maximum 
branching ratios expected. It is not enough to compute the FCNC branching ratios in, say 
the MSSM, and then evaluate them in some favorable region of the parameter space, for 
one has to preserve at the same time the stringent bounds on other observables in which 
the same physics can be applied, like the aforementioned low-energy b — > s^y decay. This 
kind of correlated study was done very carefully in Ref. [26] for the specific Higgs boson 
FCNC decays into bottom quarks within the MSSM, h -> bs (h = h°,H°,A°). In this 
paper we extend the latter work by computing also the top quark Higgs boson FCNC decay 
modes of the heavy MSSM Higgs bosons, h —> tc (h = H°, A ) under the same restrictions 
(recall that h° cannot participate in this decay because m^o < mt in the MSSM [22]). 
Furthermore, in this work we carry out an additional step absolutely necessary to make 
contact with experiment, namely we combine the FCNC decay branching ratios of the 
MSSM Higgs bosons (into both top and bottom quarks) with their MSSM production 
cross-sections in order to estimate the maximum number of FCNC events expected at 
LHC energies and luminosities. Only in this way one can assess in a practical way the 
probability of detecting such processes at the LHC. While this computation is already in 
the literature for the general 2HDM [27], to the best of our knowledge the corresponding 
calculation in the MSSM case is not available. In this paper we perform this calculation 
and compare the FCNC results obtained for the MSSM and 2HDM scenarios. 

The relevant observable quantity on which we shall focus hereafter is the cross-section 
for the production of electrically neutral pairs of heavy quarks of different flavors at the 
LHC, whose origin stems from the FCNC decays of the neutral Higgs bosons of the MSSM, 
h = h°, H°, A [22,39]. Thus, we aim at the quantity 

a(pp —>■ h — > q q') = a(pp — > hX)B(h — > qq') = cx(pp — > hX)- 



T(h -> X) 

- ( r Y ,T(h^qq' + qq') ,_ 

= a(pp hX) — — (qq = bs or tc) . (1) 

Here Y{h — > X) is the - consistently computed - total width in each case. In order to 
asses the possibility to measure these processes at the LHC, we have performed a scan 
of the MSSM parameter space to find the maximum possible value of the production 
rates (JIJ under study. The computation of the combined production rate is necessary, 
since the correlations among the different factors are important. For example, in Ref. [26] 
it was shown that the maximum branching ratio for the lightest MSSM Higgs boson, 
B(h° — > bs), is obtained in the regions of the parameter space where the coupling h° bb 
is strongly suppressed by quantum effects. On the other hand, the associated production 
a(pp — > h°bb) is one of the leading processes for the production of the lightest MSSM Higgs 
boson. It is clear then that, in the regions where B(h° — * bs) is largest, a(pp — > h°bb) will 
be suppressed. Therefore, the maximum FCNC production rate at the LHC can only be 
obtained by the combined analysis of the two relevant factors in (0) (viz. the branching 
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ratio and the Higgs boson production cross-section). We will see that the effects from 
each factor are different in different regions of the parameter space. Moreover, the re- 
alistic production FCNC rates (JTJ) in the MSSM parameter space can be obtained only 
by including the restrictions imposed by the simultaneous analysis of the branching ratio 
of the low-energy process b — > sj, whose range of values is severely limited by experi- 
ment [40-45]. As in Ref. [26], in this paper we limit ourselves to supersymmetric FCNC 
interactions mediated by the strongly interacting sector of the MSSM, i.e. the SUSY-QCD 
flavor-violating interactions induced by the gluinos. The corresponding analysis for the 
electroweak supersymmetric FCNC effects requires a lengthy separate presentation, and 
will be reported elsewhere [46]. 

The paper is organized as follows. In Section 2 we describe the general setting for 
our numerical analysis. In Section 3 we present the LHC production rates of Higgs boson 
decaying into bottom quarks through supersymmetric FCNC interactions. In Section 4 
we present the corresponding FCNC rates for the top quark channel. Finally, in Section 
5 we compare the MSSM results with the 2HDM results, and deliver our conclusions. 

2 General setting for the numerical analysis 

To compute the full one-loop value of the FCNC cross-sections cr(pp —^h^qq') for 
the three MSSM Higgs bosons (h = h°, H°, A ) we shall closely follow the notation and 
methods of Refs. [26-29]. We refer the reader to these references for the technical details. 
In particular, a thorough exposition of the relevant interaction Lagrangians and similar 
set of Feynman diagrams for the FCNC interactions is provided in [29]. See also [22,39] for 
basic definitions in the MSSM framework and [11] for detailed computational techniques 
and further illustration of the supersymmetric enhancement effects in other relevant Higgs 
boson processes. There is no need to go through these lengthy details here and in what fol- 
lows we shall limit ourselves to present the final results of our numerical analysis together 
with a detailed discussion, interpretation and phenomenological application. 

We have performed the calculations with the help of the numeric programs HIGLU, 
PPHTT [47-50] and Loop Tools [51-53]. The calculation must obviously be finite without 
renormalization, and indeed the cancellation of UV divergences using either dimensional 
regularization or dimensional reduction - the two methods giving the same results here - 
in the total amplitudes was verified explicitly. In the following we will detail the approx- 
imations used in our computation: 

• We include the full one-loop SUSY-QCD contributions to the FCNC partial decay 
widths T(h — > qq') in the observable (|T|). 

• We assume that FCNC mixing terms appear only in the LH-chiral sector of the 6x6 
squark mixing matrix. Therefore, this matrix has only non-diagonal blocks in the 
LH-LH sector. This is the most natural assumption from the theoretical point of 
view [54], and, moreover, it was proven in Ref. [29,32] that the presence of FCNC 
terms in the RH-chiral sector would enhance the partial widths by a factor two at 
most - not an order of magnitude. 
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• The Higgs sector parameters (masses and CP-even mixing angle a) have been 
treated using the leading m t and mj tan (3 approximation to the one-loop result [55- 
58]. For comparison, we also perform the analysis using the tree-level approxima- 
tion. 

• The Higgs bosons total decay widths F(h — > X) are computed at leading order, 
including all the relevant channels: T(h — > ff, ZZ, W + W~, gg). The off-shell decays 
T(h — > ZZ*,W ± W T *) have also been included. This is necessary to consistently 
compute the total decay width of T(h° —>■ X) in regions of the parameter space 
where the maximization of the cross-section (£[]) is obtained at the expense of greatly 
diminishing the partial decay widths of the two-body process h° — > bb (due to 
dramatic quantum effects that may reduce the CP-even mixing angle a to small 
values [59]). The one-loop decay rate T(h — > gg) has been taken from [49] and the 
off-shell decay partial widths have been recomputed explicitly and found perfect 
agreement with the old literature on the subject [60]. 

• The MSSM Higgs boson production cross-sections at the LHC have been computed 
using the programs HIGLU 2. 101 and PPHTT 1 . 1 [47-50]. These programs include 
the following channels: gluon-gluon fusion u(pp(gg) — > h), associated production 
with top-quarks a{pp — > hti) and associated production with bottom-quarks a(pp — > 
hbb). In order to have a consistent description, we have used the leading order 
approximation for all channels. The QCD renormalization scale is set to the default 
values for each program, namely (1q = rrih for HIGLU and fio = (rrih + 2Mq)/2 for 
PPHTT. We have used the set of CTEQ4L Parton Distribution Functions [61]. 

Running quark masses (m q (Q)) and strong coupling constants (a s (Q)) are used through- 
out, with the renormalization scale set to the decaying Higgs boson mass in the decay 
processes. More details are given below, as necessary. 

Using this setup, we have performed a maximization of the FCNC cross-section, 
Eq. (^), in the MSSM parameter space with the following restrictions on the parame- 
ters: 
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Here m g - are the LH-chiral and RH-chiral squark soft-SUSY-breaking mass parameters, 
rriq L R , common for the three generations; M g - are the physical masses of the squarks, and 
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rrih is the mass of the decaying Higgs boson h = h°, H°, A . These masses are fixed at the 
tree-level by the values of (tan/3, m^o) and the SM gauge boson masses and couplings [39]. 
Due to the structure of the Yukawa couplings in the MSSM, the value of tan/? is fixed 
at a high (small) value for the bottom (top) quark channel as indicated. The parameter 
rriAo (the mass of the CP-odd Higgs boson) is assumed to vary in the range indicated 
in (0) • At one loop these masses receive corrections from the various SUSY fields, and 
therefore depend on the values of the remaining parameters in Eq. @. The characteristic 
SUSY mass scale Msusy defines the typical mass of the squark and gluino masses 2 . The 
rest of the parameters of the squark sector are determined by this setup. For instance, by 
SU(2) gauge invariance we have m fiL = m^. Following the same notation as in [29], the 
parameter #23 represents the mixing between the second and third generation of LH-chiral 
squarks. Let us recall its definition: 

2 

5 23 = — , (3) 

m~ b m~ SL 

m~ being the non-diagonal term in the squark mass matrix squared mixing the second 
and third generation of LH-chiral squarks - and an equivalent definition for the up-type 
quarks. 3 The parameter <5 23 is a fundamental quantity in our analysis as it determines 
the strength of the tree-level FCNC interactions induced by the supersymmetric strong 
interactions, which are then transferred to the loop diagrams of the Higgs boson FCNC 
decays in Eq. (JTJ). The last two restrictions in Eq. (J2J) ensure that the (heavy) Higgs 
boson decay channels into a pair of squarks are kinematically forbidden. We have checked 
explicitly for some of the heavy Higgs boson channels (h° can never do it in practice) 
that we obtain the same results if we remove these conditions and include the partial 
widths F(h — > qq*) in the denominator of (JTJ). Strictly speaking this condition could be 
implemented, in the case of the H° boson, by just requiring m H o < 2 Mg., but we have 
made it stronger by including an additive term. This term is arbitrary (provided it is not 
very small) and acts as a buffer, namely it impedes that by an appropriate choice of the 
squark masses we can approach arbitrarily close the threshold from above, and therefore 
avoids artificial enhancement effects in our loop calculations (see below). Similarly, for the 
CP-odd Higgs boson, the condition expressed in (J2J) ensures that we avoid a similar kind 
of enhancement. In this case, however, the condition is a bit different because the A qq* 
vertex can only exist with squarks of different chirality types (A° q~ L q~ R *) or, equivalently, 
with different mass eigenstates {A° q^j*). We have used fixed values for the soft-SUSY- 
breaking trilinear couplings A t and Af, for the bs and tc channels respectively. Our results 
are essentially independent of these values and their signs. 

The task of scanning the MSSM parameter space in order to maximize a(pp — > h — > 
qq') for the various Higgs bosons is quite demanding and highly CPU-time consuming, 

2 Our programs are able to deal with completely arbitrary masses for each squark, but we are forced 
to make some simplifications in order to provide a reasonable analysis within a manageable total CPU 
time, see below. 

3 Recall that the 5ij parameters in the up-sector are related to the corresponding parameters in the 
down-sector by the Cabibbo-Kobayashi-Maskawa matrix, see e.g. [62,63]. 
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even under the conditions imposed in Eq. As stated in the introduction, our code 
includes also the restrictions on the MSSM parameter space due to the experimental 
constraint on B(b — > s-y), and therefore contains the full one-loop SUSY-QCD amplitude 
for b — > S7 constructed from the FCNC interactions induced by the gluinos. The scan 
was carried out with the help of two entirely different methods. In the first method 
we used a systematic procedure based on dividing the parameter subspace (J2J) into a 
lattice which we filled with points distributed in a completely homogeneous way. The 
second is a Monte-Carlo based method, first proposed in [64]. We have adapted the well- 
known Vegas integration program [65] to generate a sufficient number of "interesting" 
points in our parameter subspace. The total number of points used in this case was far 
smaller than in the first method. Obviously the lattice procedure gives more accurate 
results by increasing arbitrarily the total number of points, but the CPU time becomes 
prohibitively long for the whole analysis. This is so even after factoring out in a suitable 
way the phase-space integrals of the Higgs boson production processes, so that these 
integrals are computed only once for every fixed Higgs boson mass and for all the MSSM 
points of our scan in the parameter subspace (J2J). The second method is comparatively 
much faster, but it still involves a quite respectable amount of CPU time for the whole 
analysis. We found that the partial results obtained by the two methods are compatible 
at the level of 10 — 20%. For the study of our FCNC processes we consider that this level 
of accuracy should be acceptable, and for this reason all of the plots that we present in 
this work have been finally computed with the Vegas-based procedure. This also explains 
the wiggling appearance observed in the profiles of the curves presented in Sections 3-4. 
For any given abscissa point in each one of these curves, the corresponding value on the 
vertical (ordinate) axis is somewhere within a band whose width lies around 10 — 20% of 
the central value. 

A few words on the effects of the B(b —>■ sj) constraint in our analysis are now in order. 
The SUSY-QCD contribution to B(b — > s-y) can be quite large, in fact as large as the SM 
one, and with any sign. This raises the possibility of "fine-tuning" between the two type 
of contributions in certain (narrow) regions of the parameter space. As a consequence we 
could highly optimize our FCNC rates in these regions without being in conflict with the 
experimentally measured B(b —>■ 57) band. We have checked that in these regions the 
number of FCNC events can be artificially augmented by one or two orders of magnitude. 
Our scanning procedure indeed finds automatically these fine-tuning domains. However, 
we have systematically avoided them in the presentation of our analysis (for more details 
cf. [26]). In all of our plots, therefore, we show the results obtained for the non-fine-tuned 
case only. We adopt B(b — > sj) = (2.1 — 4.5) x 10~ 4 as the experimentally allowed range 
within three standard deviations [24]. 

3 Analysis of the bottom- strange channel 

The main result of the numerical scan for the bottom channel is shown in Fig. ^ To 
be specific: Fig. [T^i shows the maximum values of the production cross-sections a(pp — > 
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Figure 1: Maximum SUSY-QCD contributions to a{pp —* /i — > bs), Eq. ((H), as a function 
of (a) m^o (at fixed tan/?) and (b) tan/3 (at fixed 77140). In each plot the left-vertical axis 
provides the cross-section in pb and the right-vertical axis tracks the number of events 
per lOOfb -1 of integrated luminosity. 



h — > bs) for the three MSSM Higgs bosons h = h°, H°, A at the LHC, as a function of 
m^o; Fig. displays the cross-section as a function of tan/?. In this plot we indicate the 
value of cr(pp h — > bs) (in pb) in the left- vertical axis, and at the same time we track 
number of FCNC events (per lOOfb -1 of integrated luminosity) on the right vertical axis. 
Looking at Fig. ^ one can see immediately that at large tan/?: i) the maximum number 
of events is remarkably high (10 6 events!) for a FCNC process; ii) there is a sustained 
region in the h° channel, comprising the interval 300 GeV < m^o < 900 GeV, with a flat 
value of 5 x 10 3 events; hi) the chosen value of tan/? = 50 is not critical for H° and A 
as long it is larger than 10. In Fig. we see that the dependence on tan /? is essentially 
the same for H° and A , but for h° it is quite different: in the region 10 < tan/? < 30 the 
cross-section remains below 10 _2 pb, but for tan/? > 30 it starts climbing fast up to 0.3 pb 
at tan/? = 50. The number of events here reaches a few times 10 4 for all channels (for 
fixed mpfi = 200 GeV) 4 . For further reference, in Tabled we show the numerical values 
of a(pp — > /i — > 6s) together with the parameters which maximize the production for 
tan/? = 50 and m^o = 200 GeV. We include the value of B(h — > bs) at the maximization 
point of the FCNC cross-section. We notice that at this point the lightest Higgs boson h° 
is the one having the smallest branching ratio. This is in contrast to the situation when 
one maximizes the branching ratios independently of the number of events [26]. 

In Fig. E^, we show the effect on a(pp — > h° — ► bs), and on the total decay width 
T(h° — > X), of using the Higgs boson sector at the tree- level or at one- loop in our com- 
putation. It is well-known that the h° couplings to quarks are particularly sensitive to 

4 As already advertised, in reading the plots and tables in this work, one must keep in mind that they 
are the result of a Monte-Carlo sampling near the region of the maximal values. 
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Table 1: Maximum value of a(pp —* h bs) (and of the number of bs events per 100 fb -1 ) 
in the LHC, for m^o = 200 GeV and tan/5 = 50. Shown are also the corresponding values 
of the relevant branching ratio B(h — > bs) and of the total width of the Higgs bosons, 
together with the values of the SUSY parameters. The last row includes B{b — > sj). 

this issue, and for this reason we focus on the lightest CP-even Higgs boson for these 
considerations. We can appreciate the correlations among the different factors that enter 
the production rate (JIJ. The plotted values for a(pp — > h° — > bs) and T(h° — > X) in 
Fig. correspond precisely to the parameters that maximize (0) at the tree-level or at 
one loop in each case. Fig. shows a comparison of the various /i°-production mecha- 
nisms for the values that maximize a(pp — > h° — > 6 s) at one-loop. Remarkably, the effect 
of the radiative corrections in the Higgs sector amounts to an enhancement of our maxi- 
mal FCNC rates of up to three orders of magnitude. As we mentioned in the introduction 
(see also Ref. [26]) the maximum of B(h° — > bs) is attained under the conditions of the 
so-called "small a e g scenario" [59,66], where the two-body decay h° — > bb is strongly 
suppressed due to a corresponding suppression of the h° bb coupling. Since T(h° — > bb) 
usually dominates the total width T(h° — > X), the latter also diminishes drastically (at 
the level of the partial width of a h° three-body decay, as mentioned above). In this 
scenario the production cross-section a{pp — > h°bb) is also suppressed, so the final result 
is a compromise between the suppression of T(h° — > bb) and the possible enhancement 
(or at least sustenance) of a(pp — > h° X) by other mechanisms other than the associated 
production with bottom quarks (like the mechanism of gluon-gluon fusion, see Fig. Eb). 

Indeed, for m^o < 300 GeV we can check in Fig. that the most relevant factor for 
maximizing the FCNC cross-section is the enhancement of the h° production channel in 
association with bottom quarks, cr(pp — > h°bb). This channel operates through the bb- 
fusion vertex bb — > h° and is highly enhanced at large tan j3. In the region m^o < 300 GeV 
stays as the dominant mechanism for h° production, although one can see that the alter- 
nate gg-fusion mechanism remains all the way non-negligible. The corresponding effect 
on our FCNC cross-section ([TJ is nevertheless not obvious because this same parameter 
choice does also maximize the total width of h°, mainly through the enhancement of 
T(h° —>■ bb). There is a delicate interplay of various factors here. In particular, in the 
region m^o < 300 GeV the maximized partial width of the FCNC process h° — > bs (which 
is a function of all the parameters in (J2J)) is larger than in the region m^o > 300 GeV, 



9 



10 
10 



1 - 



> 

O 10 

u 

io" : 



10" 



10" 




r(h -> x) 

a(pp -> h° -> b s) 
r(h°-> X) tree-Higgs 



' ' ' a(pp -> h -> b s) tree-Higgs 



tanP = 50 




10 



10 



10" 



,n 4 'S 

10 g 



10 



10 



10 



200 



400 



600 



800 



100 



m A o [GeV] 



10 f 
10 3 

io 2 

10 

1 

io 1 

io" 2 
io" 3 



10" 



I 1 I 
— o(qq,gg->h bb) 
■ - ■ a(gg -> h ) 
-- a(qq,gg -> h° tt) 



tanP = 50 




200 



400 



600 



800 



1000 



m A o [GeV] 



(b) 



Figure 2: (a) h° production cross-section and decay width as a function of m^o with 
the Higgs mass relations at tree-level and at one-loop, (b) Different contributions to the 
h° production cross-section as a function of m^o corresponding to the maximization of 
a{pp — > h° — > bs) using the one-loop Higgs mass relations. 



but at the same time the total width becomes smaller in the latter region. Overall, the 
result is that a(pp — > h° — > bs) is larger in the former m^o range than in the latter. 
Furthermore, when we cross ahead the limit m^o — 300 GeV the dominant /i°-production 
channel changes turn: the associated /z°-production with bottom quarks falls abruptly 
down (see explanations below) and the gluon-gluon fusion mechanism takes over, so that 
in this range the /i°-production cross-section becomes completely dominated by gg — > h°. 
We note that while this mechanism is not so efficient at its maximum as the associated 
production, it has a virtue: it is non-suppressed in the entire range of m^o. As a re- 
sult, in the region above m^o > 300 GeV a small value of the h° bb coupling enhances 
B(h° — > bs) while it does not dramatically suppress the total cross-section a(pp — > hX). 
For m^o > 300 GeV our sampling procedure finds the maximum by selecting the points 
in parameter space corresponding to the small a e g scenario, where B(h° —>■ bs) takes 
the highest values and T(h° — > X) and a(pp —>■ h°bb) are strongly suppressed - the latter 
staying below the associated Higgs boson production with top quarks! In this way the net 
Higgs boson cross-section a{pp —* h°) is not drastically reduced thanks to the sustenance 
provided by the gg-fusion channel. Most of the loop contributions to it come from the top 
quark because the bottom quark contribution is suppressed and the squarks are rather 
heavy. One can clearly see this sustenance feature in Fig. |21 in the form of a long cross- 
section plateau up to around m^o = 1 TeV, beyond which the small a e g scenario cannot 
be maintained and T(h° — > X) starts increasing and at the same time a(pp — > h° — > bs) 
starts decreasing. It is remarkable that this behavior is only feasible thanks to the large 
radiative corrections in the MSSM Higgs sector. When we, instead, perform the computa- 
tion using the tree-level relations for the Higgs sector, the small a e g scenario is obviously 
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not possible and the enhancements/suppressions of a(pp h°bb)/T(h — > X) cannot 
take place. As a result the FCNC rate is some 3 orders of magnitude smaller than in the 
previous case. 

Next, we turn our attention to the FCNC mixing parameter $23 in Eq. (jSJ). The value 
of S23 at the cross- sect ion maximum is not necessarily the maximum allowed value of #23 
in (J2J). This is because it is a conditioned maximum, namely a maximum obtained under 
the restrictions imposed by b — » 37, as illustrated in Fig. El For further reference in our 
discussion, and to better grasp some qualitative features of our results, let us write the 
general form of the SUSY-QCD contribution to b — > 57. If we emphasize only the relevant 
supersymmetric terms under consideration (obviating the powers of the gauge couplings 
and other factors) we have 



B(b — ► S7) ~ 5. 



2 mKfjL — Abtan/3) 2 

23 77i 

JW SUSY 



(4) 



Fig.|Hk shows the maximum of a(pp — > h — > bs) as a function of ^3 for a fixed value of the 
CP-odd Higgs boson mass m^o = 200 GeV, whereas Fig. shows the computed value of 
B(b —> S7) corresponding to the parameter space points where each maximum is attained. 
At small #23 the SUSY-QCD contribution to B(b — > 37) is negligible, and the experimental 
restriction B(b —>■ 57) = (2.1 — 4.5) x 10~ 4 does not place constraints on the other MSSM 
parameters (J2J); in other words, in this region the dependence is a{pp — » h — » bs) oc (S23) 2 
- the naively expected one. Here 623 < 10 -1 ' 5 ~ 0.03 and the computed B(b — > 57) value 
lies well within the experimental limit. For larger #23, B(b 57) can be saturated at its 
uppermost experimentally allowed limit, and the rest of the parameters in (J2J) must change 
accordingly in order not to cross that limit. This can be appreciated in Fig. 0] where we 
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show the range of values taken by the physical down-type squark masses 5 (Fig. HJi) and 
the lagrangian parameters (Fig. |3t>) from Eq. (J2J) that provide the maximum values of 
a(pp — > h° — > 6 s) in Fig. El In the small region (<5 23 < 10 _L5 ~ 0.03) the parameters 
and masses that maximize a(pp —>■ h — > bs) are constant - except for the statistical 
noise unavoidable in a Monte-Carlo procedure. Note also that there are two possible 
values for the parameters \i and Aj,, due to the fact that (the leading contribution of) 
a(pp — > h — > bs) is independent of the sign of these parameters and the Monte-Carlo 
procedure picks either sign for each point with equal probability. At — 10~ L5 the 
value of B{b — > sj) becomes saturated and a(pp —> h° bs) reaches its maximum (cf. 
Fig. Eb); however #23 can keep growing, yet without overshooting the B(b — > 57) limits, 
because the increasing value of 823 is compensated by the growing squark masses (cf. 
Fig. Hk). But this is not all that simple, the higher range of ^3 can be further divided in 
two more segments where different dynamical features occur. In the first range, namely 
lQ-i-5 < £ 23 < io -0 - 75 , the heavy Higgs boson channels keep on increasing their FCNC 
rates, but not so the lightest Higgs boson channel a{pp — > h° — > bs), the reason being 
that for higher squark masses we reach the region where the small a e g scenario is feasible 
and hence the h° couplings become weakened. The relevant terms of the cross-section can 
roughly be written as follows (see Ref. [26]): 

a(pp ^bs) ~ a(pp -> h°) x ^ 3 cos 2 (/3 - a cS ) mj fi 2 /Ml VSY , (5) 

so that for large tan/5 and small a e s it becomes reduced. In the second high range of 
623, i.e. for 623 > 10~°' 75 ~ 0.18, the SUSY mass parameter Msusy has already reached 

5 There are six different down-type squarks, but four of them are nearly degenerate in mass in our 
approximation. For the 6s-channel, the down squarks are so heavy that the conditions required in the 
last two rows of @ are automatically satisfied by them in practically all the allowed range for m^o . 
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Figure 5: Maximum SUSY-QCD contributions to a(pp — > h — > bs), Eq. as a function 
of Msusy- 



its allowed maximum value specified in (j2J), therefore other parameters have to change to 
compensate for the larger 623- This is confirmed in Fig. 0Jd, where for 623 > 10 -0 ' 75 the 
absolute value of \i decreases to preserve the B{b — > 57) upper bound. Correspondingly, in 
this region a(pp —>■ h° — > bs) further falls down, as it is patent in FigEk- This additional 
feature can also be understood from the approximate expression of the cross-section given 
above. At the same time the FCNC rates for h = H°, A keep further growing, but at a 
much lower pace. This is because their (approximate) contribution goes like 

o{pp - (H°, A ) - bs) ~ a{pp - H°, A ) x 6 2 23 (sin 2 (/5 - a eS ), 1) /^7M S 4 USY , (6) 

similar to the h° case but with angular dependences on a eff and j3 which are non- 
suppressing in this region (see Ref. [26]) 6 . Again, any further increase of ^3 is now 
partially cancelled by the b — *> S7 constraint, which demands smaller values of /x. This 
explains the stabilization of the FCNC rates of the H° and A channels in the highest 
#23 range (cf. Fig. EK). The profile of the squark mass curves in Fig. 0^ implies a mixing 
mass matrix with constant diagonal terms (with value Msusy) and growing mixing terms 

(523). 

We finish our analysis of a(pp — > h — > bs) by looking at its behavior as a function of 
the SUSY mass scale Msusy > y i z - the overall scale for the squark and gluino masses - 
cf. Eq. (J2J). Fig. shows the maximum of a(pp —>■ h —>■ bs) as a function of Msusy for 
fixed m^o = 200 GeV and tan f3 = 50. The interpretation of this figure follows closely the 
results of the previous ones. At small values of Msusy the potentially large contribution 
to B{b — > sj) has to be compensated - see Eq. (jlj) - by small values of £23 and/or 
resulting in a (relatively) small value of a(pp — > h — > bs). As Msusy grows, #23 and \n\ 
can take larger values without disturbing the restrictions from B(b — > 57). The leading 
contribution to our FCNC cross-sections is actually independent of the overall SUSY mass 

6 In Eq.© above we have corrected a typing mistake present in Eq. (3.5) of Ref. [26]. 
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scale MsusY; because for all Higgs boson channels we have found the general behavior 
(leaving aside other terms mentioned above) 

<j(pp^h^bs)~<j(pp^h)x 5 2 23 m\ fi 2 /Ml VSY (h = h°,H°,A°). (7) 

The last factor effectively behaves as 5| 3 /i 2 /M| USY , and grows as 5| 3 for increasing M SUSY 
at fixed ratio /z/M S usy- Under the same conditions B(b — > 57) causes no problem because 
it is additionally suppressed by ^ft/M| USY - cf. Eq. (jlj). Therefore, we are led to a sort of 
"non-decoupling behavior" of the FCNC rates with increasing Msusy- In other words, we 
find that for the heavy neutral Higgs bosons (H°, A ) the interesting region is (contrary to 
naive expectations) the high Msusy range! The lightest Higgs boson (h°) channel shows 
a similar overall behavior, but it presents additional features because it is more tied to 
the evolution of the CP-even mixing angle a. The most interesting region for this channel 
is (cf. Fig. EJ) the central squark mass scale Msusy ~ 600 — 800 GeV, where the small a e fr 
scenario can take place. 

From the combined analysis of Figs. WS^ we arrive at the following conclusions con- 
cerning the bs final state: 

• A significant event rate of FCNC Higgs boson decays a(pp — > h — > bs) is expected 
at the LHC, even after taking into account the limits on B(b — ► sj); 

• Lightest Higgs boson case, h°: 

— For m^o < 300 GeV the rate a(pp — > h° — > bs) decreases with m^o but it is the 
largest in this interval, being produced by the combination of a large production 
cross-section a{jpp — > h°bb) and a moderate B(h° — > bs). It amounts to a 
number of events between ~ 5 x 10 3 and ~ 12 x 10 5 for every lOOfb -1 of 
integrated luminosity at the LHC; 

— For 300 GeV < m^<> < 850 GeV we expect a maximum of ~ 6 x 10 3 events/100 fb 
in the small a e ff scenario, provided by a large B(h° — » bs) and cr(pp(gg) — > h°) 
as the dominant production cross-section; 

— For rrijsfi > 850 GeV the number of events starts to decrease slowly; 

— In all cases, this maximum is attained for a large value of tan/5 ~ 50, a 
moderate value of the SUSY mass scale (Msusy ~ 600 — 800 GeV) and a low 
value of 5 23 ~ 10" 1 - 3 ~ 0.05; 

• Heavy Higgs bosons, H°,A°: 

— Although not shown in our plots, we have checked that their production rate 
a{pp — > H° A ) decreases fast with the Higgs boson mass (due to the decreasing 
of the production cross-section). We find a maximum FCNC rate of ~ 5 x 10 4 
events for m A o ~ 200 GeV, and 20 events for m^o ~ 1 TeV. 

— The maximum is produced at i) large tan/3 > 30, ii) at the highest allowed 
values of the SUSY mass scale, M S usy ~ 1 TeV, and iii) at a relatively large 
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Figure 6: Maximum SUSY-QCD contributions to a(pp — > h — > tc), Eq. as a function 
of (a) m^o (at fixed tan/3) and (b) tan/3 (at fixed m^o). 



value of the FCNC mixing parameter, #23 ~ 10 ~ 0.18, but not at the 
largest allowed value. The small a e s scenario plays no role in the heavy Higgs 
boson channels. 

Altogether one should expect a total maximum of some 120,000 events/100 fb" 1 . 



4 Analysis of the top-charm channel 

The results of the numerical scan for this channel are similar to the bs channel, so we 
will focus mainly on the differences. Fig. |HK shows the maximum value of the production 
cross-section a{pp — ► h — > tc), Eq. (JTJ), under study as a function of m^o; Fig. [Ht 1 displays 
the cross-section as a function of tan /3. Obviously the lightest Higgs boson (h°) channel 
does not appear in these plots, since in the MSSM this boson is always lighter than the 
top quark. Looking at Fig. |U] one can see immediately the following: i) the dominant 
channel in this case is the heavy scalar Higgs boson, H°; ii) it varies between 1 and 300 
events/100 fb" 1 ; iii) the tan/3 value is critical with preference for low values. In Tabled 
we show the numerical values of a(pp — > h —>■ tc) together with the parameters which 
maximize the production for tan /3 — 5 and m^o = 300 GeV. We have included the value 
of B(h —* tc) at the maximization point of the FCNC cross-section. It is remarkable that 
for the heavy CP-even Higgs boson one can reach B(H° — > tc) ~ 10~ 3 compatible with 
the b — > 57 constraint. 

Let us remark that for the heavy Higgs boson channels the features of the small a e s 
scenario play no significant role because the partial widths into bb are proportional either 
to cos 2 a e fj (in the H° case) or to sin 2 /3 (A case). Moreover, in the low tan/3 > 1 region 
(the relevant allowed one for the tc channel) the small a e g scenario does not even have a 
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Table 2: Maximum value of a(pp — > /i — > £ c) (and of the number of tc events per 100 fb -1 ) 
in the LHC, for m^o = 300 GeV and tan/3 = 5. Shown are also the corresponding values 
of the relevant branching ratio B(h — > tc) and of the total width of the Higgs bosons, 
together with the values of the SUSY parameters. The last row includes B{b — > sj). 

chance to take place. 

We turn now our view to the role of the B(b — > 57) restriction in Figs. [7| and |H1 Fig. [7| 
shows the maximum value of a(pp — > h — > tc) as a function of #23 for a fixed value of 
m^o = 300 GeV together with the corresponding computed value of B{b — > 57), while 
Fig. |S1 shows the values of the parameters, Eq. (J2J), that realize this maximum, together 
with the physical up- type squark masses. In this case the B(b —> sj) restriction is not 
as critical as in the bs channel, in part due to the fact that the SUSY-QCD contribution 
to B(b — > S7) is not enhanced at low tan/5. For 623 < 10 -1 ' 5 the value of B(b — > 
57) is well inside the experimental limits (Fig. there is no restriction on the rest 
parameters (Fig. Eb), the squark masses remain constant (Fig. |H^), and the maximum 
value of a{pp — > h — > tc) grows here in the naively expected way (S23) 2 (Fig. Ek)- Above 
this value (5 2 3 > 10~ L5 ) the parameters have to be adjusted to provide and acceptable 
range for B{b — ► s-y) 1 . In this region the SUSY mass M SU sy grows (Fig. IHJo), and \fi\ 
decreases, but not so fast as in the bs channel case (Fig. EJd). At the same time A t 
increases with increasing <5 2 3 > 10 -1 ' 5 . Most of the physical squark masses grow, but 
one of the up-type squarks (stop squark) can always have the minimum allowed mass - 
Eq. (J2J). In this region the observables under study grow more slowly, since their original 
5| 3 behavior is partially cancelled by the growing of Msusy with #23. In Fig. El we see, 
again, that a(pp — + h — > tc) grows with the SUSY mass scale (although in a way less 
pronounced than in the bs channel), due to the relaxation of the b — > 37 constraint for 
large Msusy- While a(pp — > h — > bs) is augmented nearly three orders of magnitude in 
the range Msusy = 200 — 1000 GeV (Fig. |5J), the tc channel undergoes only an increase 
of roughly a factor 10 in the same parameter range. 

How do the two Higgs boson channels H° and A compare as sources of tc events? 
The gg-fusion mechanism is one of the leading processes for Higgs boson production at 

7 The two lines appearing in this region for Bn{b — > sj) mean that the maximum of cr(pp — > H° — > tc) 
is attained either by the maximum or the minimum allowed value of B(b — > sj), our Monte-Carlo sampling 
procedure picks either choice with equal probability for each value of 623- 
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Figure 7: Maximum SUSY-QCD contributions to a(pp — > h — ► £c), Eq. (JTJ), as a function 
of (a) 5 2 3 and (b) value of 6 — > 57, for fixed tan/? and m^o. 
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Figure 9: Maximum SUSY-QCD contributions to a(pp — > h — > tc), Eq. (JTJ), as a function 
of Msusy- 

relatively small values of tan /3 > 1 - associated production with bb remaining still sizeable. 
Due to CP conservation, the squark contributions to gg — > A cancel out at one-loop and 
only the quark contributions remain [50]. For large squark masses the production cross- 
sections for H° and A are similar, the latter being slightly larger. We see from Fig-EI 
that, despite the similarity in production, the H° channel gives larger FCNC rates than 
the A one. The excess of FCNC events from the former can be explained mainly from 
the constraints that we have imposed from the very beginning on the squarks masses 
in relation to the Higgs boson masses (see Eq. As we have noted in section |21 we 
can have squarks of the same chirality-type in the H° qq* vertex, whereas they must 
necessarily be of opposite chirality-type in the A°qq* case. As a result, for small m^o 
the squark mass constraints expressed in (j2J) allow the FCNC cross-section maximization 
process to pick points near the saturation of the mass condition 2 Mq. > m H o + 50 GeV, 
but not of M^+Mq. > m^o + 50 GeV for {i ^ j) because only one of the up-squarks can be 
light (see Fig. |H^). This produces an enhancement of the branching ratio of H° — > tc and 
for this reason this FCNC channel dominates in the relatively small m^o region. However, 
as soon as m^o is sufficiently heavy the second mass constraint can also be satisfied and 
then the two curves in Fig. Intend to converge. 

From the combined analysis of Figs. IHK1 we conclude that in the case of the H° 
channel we expect a maximum of ~ 300 events/100 fb _1 decays into top quarks at the 
LHC. This maximum is achieved for a CP-odd Higgs boson mass of m^o ~ 300 GeV, and 
a moderately low tan/5 ~ 5. This rate can grow one order of magnitude by a lower value 
of tan/3 ~ 2, but decreases significantly with m^o. The maximum is obtained at the 
largest possible value of 5 2 3, and a moderate SUSY mass scale M SU sy ~ 600 — 800 GeV, 
but having one of the squarks light. While the number of events is significantly lower 
than the 6s-channel ones, the tc-channels offers a better opportunity for detection, due to 
the much lower background. 
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5 Discussion and conclusions 



We have carried out a systematic study of the production rate of FCNC processes at the 
LHC mediated by the decay of neutral Higgs bosons of the MSSM: a{pp — > h — > q q') (h = 
h°, H°, A ) - see Eq. ([!]). Specifically, we have concentrated on the FCNC production of 
the heavy quark pairs qq' = bs and tc, because they are the only ones that have a chance of 
being detected. We have focused on the FCNC supersymmetric effects stemming from the 
strongly interacting sector of the MSSM, namely from the gluino-mediated flavor- changing 
interactions. We have performed a maximization of the event rates in the parameter space 
under a set of conditions that can be considered "irreducible", see Eq. ©, i.e. we cannot 
further shorten this minimal set (e.g. by making additional assumptions on the relations 
among the parameters) without potentially jeopardizing the conclusions of this study. 
Even within this restricted parameter subspace the computer analysis has been rather 
demanding. The numerical scan has been performed using Monte Carlo techniques which 
we have partially cross-checked with more conventional methods. The maximization of 
the cross-sections (|TJ) has been performed by simultaneously computing the corresponding 
MSSM quantum effects on the (relatively well-measured) low-energy FCNC decay b — *■ S7 
and requiring that the experimental limits on this observable are preserved. 

To summarize our results: the total number of FCNC heavy flavor events originating 
from supersymmetric Higgs boson interactions at the LHC can be large (of order 10 6 ), but 
this does not mean that they can be easily disentangled from the underlying background 
of QCD jets where they are immersed. For example, it is well known that the simple two- 
body decay h — > bb is impossible to isolate due to the huge irreducible QCD background 
from bb dijets - a result that holds for both the SM and the MSSM [1,2]. This led a 
long time ago to complement the search with many other channels, particularly 
which has been identified as an excellent signature in the appropriate range. Similarly, 
the FCNC Higgs boson decay channels may help to complement the general Higgs boson 
search strategies, mainly because the FCNC processes should be essentially free of QCD 
background. Notwithstanding other difficulties can appear, such as misidentification of 
jets. For instance, for the 6s final states misidentification of 6-quarks as c-quarks in cs- 
production from charged currents may obscure the possibility that the 6s-events can be 
really attributed to Higgs boson FCNC decays. This also applies to the tc final states, 
where misidentification of 6-quarks as c-quarks in e.g. tb production might be a source 
of background to the tc events, although in this case the clear-cut top quark signature 
should be much more helpful (specially after an appropriate study of the distribution of 
the signal versus the background). However, to rate the actual impact of these disturbing 
effects one would need an additional study which is beyond the scope of the present paper. 

An interesting (and counter-intuitive) result of our work is that for all the Higgs boson 
channels h = h°, H°, A , the FCNC cross-section a(pp — > h — > qq') increases with growing 
SUSY mass scale Msusy- Due to this effective "non-decoupling" behavior (which is more 
pronounced for the bs channel) the FCNC rates are maximal when the overall squark 
and gluino mass scale is of order of Msusy i$ 1 TeV - with the only proviso that for 
the tc-channel a single squark should have a low mass (> 150 GeV). Moreover, we find 
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that the two types of FCNC final states (bs and tc) prefer different ranges of tan (3. The 
fes-channel is most efficient at high tan/5 > 30, whereas the tc-channel works better in 
the regime of low tan/3 < 10 (see below). As for the mixing parameter #23 (which is the 
fundamental supersymmetric FCNC parameter of our analysis, see its definition in (JSJ)) 
we remark that the maximum number of events is not always attained for the largest 
possible values of it (due to the influence of the b — > 37 constraint): in the 6s-channel 
the maximum is achieved for moderate values in the range 0.05 < #23 < 0.2, whereas 
the tc-channel prefers the maximum allowed values in our analysis (S23 < 0.8). We also 
remark that the naive expectation a(pp — > h — > qq') ~ 5| 3 does not always apply. 

The number of FCNC events originating from the two channels bs and tc is not alike, 
and it also depends on the particular Higgs boson. For the bs final states we have found 
that the optimized value of a(pp — > h — > bs) produced by our analysis is ~ 12 pb. This 
amounts to ~ 12 x 10 5 events per J Cdt = 100 fb -1 of data at the LHC. The most favorable 
Higgs boson channel is the one corresponding to the lightest MSSM Higgs boson, h°. For 
this boson there are non-trivial correlations between the two factors in Eq. (jT|. namely 
between the Higgs boson production cross-section and the FCNC branching ratio. These 
correlations permit an increase of the total number of FCNC events up to two orders of 
magnitude in certain compared to the number of events produced by the heavy 

Higgs bosons H° and A , which are essentially free of these correlations. The latter stem 
from relevant quantum effects on the parameters of the Higgs boson sector at one-loop 
precisely in the regions of our interest. 

On the other hand, the maximum value of a(pp — > h — > tc) is more moderate, to 
wit: 3 x 10" 3 pb, or ~ 300 events/100 fb . For the total integrated luminosity during 
the operative lifetime of the LHC, which amounts to some (300 — 400) fb -1 , we estimate 
that a few thousand tc events could be collected in the most optimistic conditions. This 
number is of course sensitive to many MSSM parameters, but most particularly to two: 
m^o and tan/5. The mass of the CP-odd Higgs boson should not be heavier than m^o ~ 
(400—500) GeV if one does not want to decrease the number of events below a few hundred 
per 100 fb -1 . On the other hand the number of events is very much dependent on the 
particular range of tan (3. As we have said above, the lowest possible values are preferred 
for the tc channel, but the sensitivity in this range is so high that the order of magnitude 
of a(pp — > h — > tc) may change for different (close) choices of tan/3. Throughout all the 
analysis of the tc channel we have fixed tan f3 at an intermediate value, but the maximum 
number of events per 100 AT 1 would grow from ~ 300 (for our standard choice tan (3 = 5) 
up to ~ (500,900,2000) if we would have chosen tan/3 = (4,3,2) respectively. Such 
lower values of tan/? are usually avoided in some MSSM analyses in the literature, but 
as a matter of fact there are no fully water-tight experimental bounds on tan (3 excluding 
this lower range, apart from the more incontrovertible strict lowest limit tan/? > 1. We 
recall that the lower limit on tan/? is obtained indirectly from the LEP exclusion data 
on the light Higgs boson search, and is therefore very sensitive to the inputs used in the 
computation, specially the top quark mass [67]. Unlike the difficulties in the 6s-channel, 
and in spite of the substantially smaller number of events, we deem more feasible to 
extract the tc signal at the LHC, due to the presence of the quark top, which carries a 
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highly distinguishable signature. 

At this point a comparison with our previous work [27] is in order. In that work we 
have studied in quite some detail the maximum FCNC production rates of Higgs bosons 
decaying into tc final states within the general two-Higgs-doublet model. It was found 
that the maximal branching ratio in the 2HDM takes place in the type II 2HDM (or 
2HDM II), and reads B u (h -> tc) ~ 1(T 5 , whereas in the 2HDM I it is comparatively 
negligible. After a detailed computation of the event rates (including also the particular 
restrictions of the b — > 57 process, which are different in the 2HDM case as compared to 
the MSSM) the conclusion was that several hundred tc events could be collected at the 
LHC under optimal conditions. We clearly identified which are the most relevant Higgs 
boson modes for this purpose and the domains of the 2HDM II parameter space where 
these events could originate from. To make it short, our conclusion there [27] was that the 
h° is the most gifted decay and the ideal situation occurs when tan j3 and tan a are both 
large, and also when the CP-odd state A is much heavier than the CP-even ones (h°, H°). 
Furthermore we found that in the general 2HDM the A state never gives any appreciable 
FCNC rate into tc. It is easy to see that the mode A — > bs is not favored either (unless 
tan/5 is very small and m H ± unusually light, both situations rather unappealing). 

How it compares with the MSSM case under study? To start with, we note that here 
the A channel gives essentially the same bs rate as the H° one, and that both modes 
can be quite relevant. At the same time the A rate into tc is, though not dominant, not 
negligible at all. Some few hundred events of this nature per 100 fb _1 are possible. If this 
is not enough, in the MSSM the most relevant tan f5 region for the tc final states is not the 
highest one (as in the 2HDM II case) but just the opposite: the lowest allowed one. This 
is an important difference, and one that should help to discriminate between the 2HDM 
and MSSM models in case that some tc events would be unambiguously tagged at the 
LHC. After all there are many high precision observables that are highly sensitive to the 
preferred range of tan /3, so that the favorite value of tan (3 could already be fixed from 
other experiments by the time that some FCNC events could be detected. Apart from 
the different correlation of the parameters in the non-supersymmetric and supersymmetric 
model, in the latter case the maximal event rates for the tc mode are typically one order 
of magnitude higher than the maximal event rates in the former. 

The corresponding study for the 2HDM branching ratios into bs final states was per- 
formed in [37], although in this work the cross-section and number of events were not 
computed. However, from the maximum size of the expected branching ratios compatible 
with b — > 57 (viz. B n (h -> bs) ~ 10~ 6 - 10~ 5 in the 2HDM II) it is already pretty obvious 
that the number of events can never be competitive with the supersymmetric case where 
B MSSM (h -> bs) ~ 10~ 4 - 10~ 3 [26]. As for the 2HDM I (which is insensitive to the 
b — > S7 bounds) the branching ratios can be at most of order B : {h — > bs) ~ 10~ 5 — 10~ 3 
and only so for very small values of tan/5 = 0.1 — 0.5 which are actually excluded in the 
MSSM. 

Summing up and closing: the maximum number of FCNC events in the MSSM case 
is larger than the highest expected rates both in the 2HDM I and II , the two kind of 
signatures of physics beyond the SM being perfectly distinguishable because the relevant 
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regions of the parameter space are completely different. If a sample of FCNC events of 
this kind could be collected, we should be able to ascertain which is its ultimate origin. At 
the end of the day if one single thing should be emphasized is that the FCNC event rate 
into bs or tc is so extremely tiny in the SM that if only a dozen events of this kind could be 
captured under suitable experimental conditions it would be an undeniable signature of 
new physics. From what we have seen, the odds should be heavily in favor of attributing 
it to a supersymmetric origin. 
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